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We present a protocol for the teleportation of the quantum state of a pulse of light onto the 
collective spin state of an atomic ensemble. The entangled state of light and atoms employed as 
a resource in this protocol is created by probing the collective atomic spin, Larmor precessing in 
an external magnetic field, off resonantly with a coherent pulse of light. We take here for the first 
time full account of the effects of Larmor precession and show that it gives rise to a qualitatively 
new type of multimode entangled state of light and atoms. The protocol is shown to be robust 
against the dominating sources of noise and can be implemented with an atomic ensemble at room 
temperature interacting with free space light. We also provide a scheme to perform the readout of 
the Larmor precessing spin state enabling the verification of successful teleportation as well as the 
creation of spin squeezing. 

PACS numbers: 03.67.Mn, 32.80.Qk 



I. INTRODUCTION 

Quantum teleportation - the disembodied transport of 
quantum states - has been demonstrated so far in sev- 
eral seminal experiments dealing with purely photonic 
or atomic systems. Here we propose a protocol for 
the teleportation of a coherent state carried initially by 
a pulse of light onto the collective spin state of ^ 10^^ 
atoms. This protocol - just as the recently demonstrated 
direct transfer of a quantum state of light onto atoms Q 

- is particularly relevant for long distance entanglement 
distribution, a key resource in quantum communication 
networks Q. 

Our scheme can be implemented with just coherent 
light and room-temperature atoms in a single vapor cell 
placed in a homogeneous magnetic field. Existing proto- 
cols in Quantum Information (QI) with continuous vari- 
ables of atomic ensembles and light 4] arc commonly de- 
signed for setups where no external magnetic field is ap- 
plied such that the interaction of light with atoms meets 
the Quantum non-demolition (QND) criteria [H>l^' In 
contrast, in all experiments dealing with vapor cells at 
room-temperature 0, 13 it is, for technical reasons, ab- 
solutely essential to employ magnetic fields. In experi- 
ments Q two cells with counter-rotating atomic spins 
were used to comply with both, the need for an exter- 
nal magnetic field and the one for an interaction of QND 
character. So far it was believed to be impossible to use a 
single cell in a magnetic field to implement QI protocols, 
since in this case - due to the Larmor precession - scat- 
tered light simultaneously reads out two non-commuting 
spin components such that the interaction is not of QND 
type. 

In this paper we do not only show that it is well 
possible to make use of the quantum state of light and 
atoms created in this setup but we demonstrate that 

- for the purpose of teleportation [H, Isj - it is in fact 



better to do so. As compared to the state resulting 
from the common QND interaction the application of 
an external magnetic field enhances the creation of 
correlations between atoms and light, generating more 
and qualitatively new, multimode type of entanglement. 
The results of the paper can be summarized as follows: 
(i) Larmor precession in an external magnetic field 
enhances the creation of entanglement when a collective 
atomic spin is probed with off-resonant light. The 
resulting entanglement involves multiple modes and is 
stronger as compared to what can be achieved in a 
comparable QND interaction. 

{a) This type of entangled state can be used as a 
resource in a teleportation protocol, which is a simple 
generalization of the standard protocol 's', based on 
Einstein-Podolsky-Rosen (EPR) type of entanglement. 
For the experimentally accessible parameter regime the 
teleportation fidelity is close to optimal. The protocol 
is robust against imperfections and can be implemented 
with state of the art technique. 

{iii) Homodyne detection of appropriate scattering 
modes of light leaves the atomic state in a spin squeezed 
state. The squeezing can be the same as attained from 
a comparable QND measurement of the atomic spin 
[lOl . The same scheme can be used for atomic state 
read-out of the Larmor precessing spin, necessary to 
verify successful teleportation. 



We would like to note that it was shown recently in 
T^l that the effect of a magnetic field can enhance the 
capacity of a quantum memory in the setup of two cells. 
Teleportation in the setup of a single cell without mag- 
netic field was addressed in [l3l |. 

The paper is organized as follows: The three points 
above are presented in sections m IIIII and IIVI in this 
order. Some of the details in the calculations of sections 
IIIII and Hvl are moved to appendices IbI and lO 
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II. INTERACTION 

We consider an ensemble of Nat Alkali atoms with total 
ground state angular momentum _F, placed in a constant 
magnetic field causing a Zeeman splitting of Ml and ini- 
tially prepared in a fully polarized state along x. The 
collective spin of the ensemble is then probed by an off 
resonant pulse which propagates along z and is linearly 
polarized along x. Thorough descriptions of this inter- 
action and the final state of light and atoms after the 
scatterin g ca n be found in __14, 15, 16, 17, 18] and espe- 
cially in [l^ l2(l I21I l23 | for the specific system we have 
in mind. We derive the final state here with a special 
focus on the effects of Larmor precession and light prop- 
agation in order to identify the light modes which arc 
actually populated in the scattering process. 

In appendix ^ we show that the interaction is ade- 
quately described by a Hamiltonian 
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Hat + Hu + F, 



(1) 



and where Hu is the Hamiltonian for the free radi- 
ation field. The canonical conjugate variables X, P 
describe in the Holstein-Primakoff approximation j2,l | 
transverse components of the collective angular mo- 
mentum in y and z direction respectively. They satisfy 
[X, P] ^ i and have zero mean and a normalized variance 
AX'^ — AP^ — 1/2 for the initial coherent spin state. 
In analogy to light field quadratures we will denote the 
normalized transverse spin components X, P also as spin 
quadratures. For the light field only its linearly polarized 
component along y is relevant and is described in terms 
of quadratures of spatially localized modes 0, 
x{z),p{z), which obey [x{z),p{z')] = ic6{z — z'). Before 
the interaction process, this polarization component 
is in vacuum such that initially {x{z)) = {p{z)) = 
and {x{z)x{z')) = {p{z)p{z')) = c5{z - z')/2. The 
dimensionless coupling constant is given by 
K = NphNatFai<y^ /"^AA where Nph is the over- 
all number of photons in the pulse, ai is a constant 
characterizing the ground state's vector polarizability, a 
is the scattering cross section, T the decay rate, A the 
detuning and A the effective beam cross section. 

Changing to a rotating frame with respect to Hat by 
defining Xi{t) = eKp(—iHatt)Xexp{iHatt) and evaluat- 
ing the Heisenberg equations for these operators yields 
the following Maxwell-Bloch equations 
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where dt(z) denotes the partial derivative with respect to 
t (z). These equations have a clear interpretation. Light 
noise coming from the field in quadrature with the classi- 
cal probe piles up in both, the X and P spin quadrature, 
but it alternately affects only one or the other, changing 
with a period of 1/17. Conversely atomic noise adds to 
the in phase field quadrature only and the signal comes 
alternately from the X and P spin quadrature. The out 
of phase field quadrature is conserved in the interaction. 

To solve this set of coupled equations it is convenient 
to introduce a new position variable, ^ — ct — z, to elim- 
inate the z dependence. New light quadratures defined 
by x{£^,t) = x{ct — ^,t), p{^,t) = p{ct — ^,t) also have 
a simple interpretation: ^ labels the slices of the pulse 
moving in and out of the ensemble one after the other, 
starting with ^ = and terminating at ^ = cT. The 
Maxwell equations now read 



dtpitt) = 0, 
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The solutions to equations (|2al 12 bl Ec|l are 

Xi{t) = X/(0) + -j= dTCOs(rjT)p(cr,0), (3a) 

Vt Jo 

Pi{t) = PiiO) + —= / dr sin(f7T)p(cr, 0) , (3b) 
VT Ja 



Pitt)^pii,0) 
and the formal solution to Ij2d|l is 
x{tt)^x{tO) + 



(3c) 



(3d) 



[cosine /c) Pi (a c) - sm{ne/c)XiiC/c)] 



As mentioned before, both atomic spin quadratures 
are affected by light but, as is evident from the solutions 
for X{t), P{t), they receive contributions from different 
and, in fact, orthogonal projections of the out-of-phase 
field. As we will show in the following, the corresponding 
projections of the in-phase field carry in turn the signal 
of atomic quadratures after the interaction. It is there- 
fore convenient to explicitly introduce operators for these 
modes .21j . We define a cosine component before the in- 
teraction 
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dr cos(r2r)x(cr, 0) 



(4a) 



(4b) 



and a sine component p^", x™ with cos(r2T) replaced by 
sin(riT). In frequency space these modes consist of spec- 
tral components at sidebands Wc ± and are closely re- 
lated to the sideband modulation modes introduced in 
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[26] for the description of two photon processes. It is eas- 
ily checked that these modes are asymptotically canon- 
ical, [x)^,p'^] = [a;^",pf ] = i[l + 0{n~^)] ~ i, and inde- 
pendent, — ©(tiq ^) ~ 0, if we assume no ^ 1 
for no — QT, the pulse length measured in periods of 
Larmor precession. 

In terms of these modes the atomic state after the in- 
teraction X™* = Xi{T), P™* = Pi{T) is given by 
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The final 



state of 

out ^out 



scribed by x-^^^,p^(^^^ 



cosine (sine) modes is de- 
defined by equations Q with 
a;(cr, 0), p(cT, 0) replaced by x(ct,T), p[cT,T) respec- 
tively. Since the out-of-phase field is conserved we have 
trivially 
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FIG. 1: Von Neumann Entropy of the reduced state of atoms 
versus coupling strength kappa for the state of equation |3 
(full line) and for the state generated without magnetic field 
in a pure QND interaction (dashed line) with the same cou- 
pling strength. Application of a magnetic field significantly 
enhances the amount of light-atom entanglement. 



(5b) 



Deriving the corresponding expressions for the cosine and 
sine components of the field in phase, a;™*, raises 

some difficulties connected to the back action of light 
onto itself. This effect can be understood by noting that 
a slice £, of the pulse receives a signal of atoms at a time 
^/c [see equation (j3d|l ] which, regarding equations 
13 b|) . in turn carry already the integrated signal of all 
slices up to ^. Thus, mediated by the atoms, light acts 
back on itself. The technicalities in the treatment of this 
effect are given in appendix^where we identify relevant 
"back action modes", x ^ i, p ^ x ^ i, p ^ in terms of 
which one can express 
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The last two terms in both lines represent the effect of 
back action, part of which involves the already defined 
cosine and sine components of the field in quadrature. 
The remaining part is subsumed in the back action modes 
which are again canonical and independent from all other 
modes. 

Equations Q describe the final state of atoms and the 
relevant part of scattered light after the pulse has passed 
the atomic ensemble and are the central result of this 
section. Treating the last terms in equations (|5cl5d() as 
noise terms, it is readily checked by means of the separa- 
bility criteria in ^2^ that this state is fully inseparable, 
i.e. it is inseparable with respect to all splittings be- 
tween the three modes. For the following teleportation 
protocol the relevant entanglement is the one between 
atoms and the two light modes. Figure ^ shows the von 
Neumann entropy £^vN of the reduced state of atoms in 
its dependence on the coupling strength n and in com- 
parison with the entanglement created without magnetic 
field in a pure QND interaction of atoms and light. The 
amount of entanglement is significantly enhanced. 



III. TELEPORTATION OF LIGHT ONTO 
ATOMS 

In this section we will show how the multimode entan- 
glement between light and atoms generated in the scat- 
tering process can be employed in a teleportation proto- 
col which is a simple generalization of the standard proto- 
col for continuous variable teleportation using EPR-type 
entangled states [3, . We first present the protocol and 
evaluate its fidelity and then analyze its performance un- 
der realistic experimental conditions. 



A. Basic protocol 

Figure El depicts the basic scheme which, as usually, 
consists of a Bell measurement and a feedback operation. 

Input The coherent state to be teleported is encoded 
in a pulse which is linearly polarized orthogonal to the 
classical driving pulse and whose carrier frequency lies at 
the upper sideband, i.e. at ujc + ^. The pulse envelope 
has to match the one of the classical pulse. As is shown 
in appendix IbI canonical operators y, q with [y, q] — i de- 
scribing this mode can conveniently be expressed in terms 
of cosine and sine modulation modes, analogous to equa- 
tions l@J, defined with respect to the carrier frequency. 
One finds 



q = --^{yc-qs)- (6) 



A coherent input amounts to having initially 
Ay^ ~ Aq^ — 1/2 and an amplitude (y), (q) with 
mean photon number Uph — ((y)^ + ('?)^)/2- 

Bell measurement This input is combined at a beam 
splitter with the classical pulse and the scattered light. 
At the ports of the beam splitter Stokes vector compo- 
nents Sy and Sz are measured by means of standard po- 
larization measurements. Given the classical pulse in x 
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Feedback 



we find by means of equations (O, © and Q 
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FIG. 2: Scheme for teleportation of light onto atoms: As de- 
scribed in section^ a classical pulse (linearly polarized along 
x) propagating along the positive z direction is scattered off 
an atomic ensemble contained in a glass cell and placed in a 
constant magnetic field B along x. Classical pulse and scat- 
tered light (linearly polarized along y) are overlapped with 
a with a coherent pulse (linearly polarized along z) at beam 
splitter BS. By means of standard polarization measurements 
Stokes vector components 5*^ and Sz are measured at one and 
the other port respectively, realizing the Bell measurement. 
The Fourier components at Larmor frequency SI of the corre- 
sponding photocurrents determine the amount of conditional 
displacement of the atomic spin which can be achieved by ap- 
plying a properly timed transverse magnetic field b{t). See 
section Hira for details. 



polarization this amounts to a homodyne detection of 
in- and out-of-phase fields of the orthogonal polarization 
component. The resulting photocurrents are numerically 
demodulated to extract the relevant sine and cosine com- 
ponents at the Larmor frequency |23|. Thus one effec- 
tively measures the commuting observables 



i=(xr + 2/c), Hs^^ixT' + ys), (7) 
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Let the respective measurement results be given by 
Xc, Xs, Qc and Qs- 

Feedback Conditioned on these results the atomic 
state is then displaced by an amount Xg — Qc in ^ and 
—Xc — Qs in P. This can be achieved by means of two 
fast radio- frequency magnetic pulses separated by a quar- 
ter of a Larmor period. In the ensemble average the final 
state of atoms is simply given by 



X*""" = X' 



This description of feedback is justified rigorously in ap- 
pendixEl Relating these expressions to input operators, 
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This is the main result of this section. 

Teleportation fidelity Taking the mean of equations 
@ with respect to the initial state all contributions due 
to input operators and back action modes vanish such 
that (X'*") = (y) and (pf*") = {q). Thus, the am- 
plitude of the coherent input light pulse is mapped on 
atomic spin quadratures as desired. In order to prove 
faithful teleportation also the variances have to be con- 
served. It is evident from that the final atomic spin 
variances will be increased as compared to the coherent 
input. These additional terms describe unwanted excess 
noise and have to be minimized by a proper choice of 
the coupling k. As a figure of merit for the telepor- 
tation protocol we use the fidelity, i.e. squared over- 
lap, of input and final state. Given that the means are 
transmitted correctly the fidelity is found to be F = 
2 [(l-|-2(AXf*")2)(l-h2(AP*i")2)]~^''^ The variances 
of the final spin quadratures are readily calculated tak- 
ing into account that all modes involved are independent 
and have initially a normalized variance of 1/2. In this 
way a theoretical limit on the achievable fidelity can be 
derived depending solely on the coupling strength k. In 
figure|3|we take advantage of the fact that the amount of 
entanglement between light and atoms is a monotonously 
increasing function of k, such that we can plot the fi- 
delity versus the entanglement. This has the advantage 
that we can compare the performance of our teleporta- 
tion protocol with the canonical one H, 0| which uses 
a two-mode squeezed state of the same entanglement as 
a resource and therefore maximizes the teleportation fi- 
delity for the given amount of entanglement. No physical 
state can achieve a higher fidelity with the same entan- 
glement. This follows from the results of where it 
was shown that two-mode squeezed states minimize the 
EPR variance (and therefore maximize the teleportation 
fidelity) for given entanglement. The theoretical fidelity 
achievable in our protocol is maximized for k ~ 1.64 cor- 
responding to P ~ .77. But also for experimentally more 
feasible values of k ~ 1 can the fidelity well exceed the 
classical limit |2^ |2^ of 1/2 and, moreover, compari- 
son with the values achievable with a two-mode squeezed 
state shows that our protocol is close to optimal. 
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FIG. 3: (a) Theoretical limit on the achievable fidelity F ver- 
sus entanglement between atoms and light measured by the 
von Neumann entropy _EvN of the reduced state of atoms. The 
grey area is unphysical. For moderate amounts of entangle- 
ment our protocol is close to optimal, (b) Coupling strength k 
versus entanglement. The dashed lines indicate the maximal 
fidelity of _F = .77 which is achieved for k = 1.64. 



B. Noise effects and Gaussian distributed input 

Under realistic conditions the teleportation fidelity will 
be degraded by noise effects like decoherence of the 
atomic spin state, light absorption and reflection losses 
and also because the coupling constant k is experimen- 
tally limited to values k ~ 1. On the other hand the 
classical fidelity bound to be beaten will be somewhat 
higher than 1/2 since the coherent input states will nec- 
essarily be drawn according to a distribution with a finite 
width in the mean photon number n. In this section we 
analyze the efficiency of the teleportation protocol un- 
der these conditions and show that it is still possible to 
surpass any classical strategy for the transmission and 
storage of coherent states of light (28^ .29] . 

During the interaction atomic polarization decays due 
to spontaneous emission and collisional relaxation. In- 
cluding a transverse decay the final state of atoms is given 
by 

= ^/T^(x- + -^p|^")-fv^/x, (lOa) 
= ^/]—p(P"^ + ^p^^) + ^fp. (lOb) 

as follows from the discussion in appendix ^ [3 is the 
atomic decay parameter and /x, fp are Langevin noise 
operators with zero mean. Their variance is experimen- 
tally found to be close to the value corresponding to a 
coherent state such that {fx) = {fp) — 1/2. 

Light absorption and reflection losses can be taken into 
account in the same way as finite detection efficiency. For 
example the statistics of measurement outcome Xg will 



not stem from the signal mode Xg alone but rather from 
the noisy mode ^/l — exg + y/efx,s where e is the photon 
loss parameter and fx.s is a Langevin noise operator of 
zero mean and variance {f^s) — 1/2- Analogous expres- 
sions have to be used for the measurements of Xc , Qs and 
Qc which will be adulterated by Langevin terms fx,c, fq,s 
and fq^c respectively. In principle each of the measure- 
ment outcomes can be fed back with an independently 
chosen gain but for symmetry reasons it is enough to 
distinguish gain coefficients gx, gq for the measurement 
outcomes of sine and cosine components of x and q re- 
spectively. Including photon loss, finite gain and atomic 
decay, as given in ((Tn|) , equations © , describing the final 
state of atoms after the feed back operation, generalize 
to 
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For non unit gains a given coherent amplitude {{y) , {q)) 
will not be perfectly teleported onto atoms and the cor- 
responding fidelity will be degraded by this mismatch 
according to 

F{{y)Aq)) 



■ exp 



2(AXfi")2][l + 2(APfi")2] 

{{y)-{X^-))' ((g)_(pfin^)2 

l + 2(AXfi")2 1-|-2(AP 



Rn\2 



If the input amplitudes are drawn ac- 
cording to a Gaussian distribution 
p{{y), {q)) — exp[— ((j/)2 + {q)^)/2n]/2TTn with mean 
photon number h the average fidelity [with respect to 
((y), (g))] is readily calculated. The exact expression in 
terms of initial operators can then be derived by means 
of equations lO, ©, and l(TT|l but is not particularly 
enlightening. In figure 01 we plot the average fidelity, 
optimized with respect to gains gx, gq, in its dependence 
on the atomic decay (3 for various values of photon loss 
e. We assume a realistic value k = 0.96 for the coupling 
constant and a mean number of photons n = 4 for the 
distribution of the coherent input. For feasible values 

(3, e 0.2 the average fidelity is still well above the 
classical bound on the fidelity ze, "2^ . This proves that 
the proposed protocol is robust against the dominating 
noise effects in this system. 

The experimental feasibility of the proposal is illus- 
trated with the following example. Consider a sam- 
ple of Nat — 10^^ Cesium atoms in a glass cell placed 
in a constant magnetic field along the a;-direction caus- 
ing a Zeeman splitting of = 350 kHz in the F = 4 
ground state multiplet. The atoms are pumped into 
mp = 4 and probed on the D2 (P = 4 ^ F' = 3, 4, 5) 
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FIG. 4: (a) Average fidelity achievable in the presence 
of atomic decay (i, reflection and light absorption losses 
e = 8%, 12%, 16%, coupling k = 0.96 and Gaussian dis- 
tributed input states with mean photon number n = A. The fi- 
delity benchmark is in this case 5/9 (dashed line), (b) Respec- 
tive optimal values for gains (solid lines) and Qq (dashed 
lines) . 



transition. The classical pulse contains an overall num- 
ber of Nph = 2.5 10^'^ photons, is detuned to the blue 
by A = 1 GHz, has a duration T — 1 ms and can have 
an effective cross section oi A ~ 6cm^ due to thermal 
motion of atoms. Under these conditions the tensor 
polarizability can be neglected (A/whfs — 10~^). Also 
riQ = = 350 justifies the use of independent scatter- 
ing modes. The coupling k ~ 1 and the depumping of 
ground state population 77 ~ 10^^ as desired. 



IV. SPIN SQUEEZING AND STATE READ-OUT 

In this section we present a scheme for reading out 
either of the atomic spin components X, P by means of 
a probe pulse interacting with the atoms in the one way 
as described in section ^ The proposed scheme allows 
one, on the one hand, to verify successful receipt of the 
coherent input subsequent to the teleportation protocol 
of section imi and, on the other hand, enables to generate 
spin squeezing if it is performed on a coherent spin state. 

It is well known [Tsl l3]| and was demonstrated exper- 
imentally [lol [lll | that the pure interaction V , as given 
in equation can be used to perform a QND measure- 
ment of either of the transverse spin components. At 
first sight this seems not to be an option in the scenario 
under consideration since the local term Hat^ account- 
ing for Larmor precession, commutes with neither of the 
spin quadratures such that the total Hamiltonian does 
not satisfy the QND criteria |^ . As we have shown in 
section ^ Larmor precession has two effects: Scattered 



light is correlated with both transverse components and 
suffers from back action mediated by the atoms. Thus, 
in order to read out a single spin component one has to 
overcome both disturbing effects. 

Our claim is that this can be achieved by a simultane- 
ous measurement of p°^i or _p°"*, p°^i 
if, respectively, X or P is to be measured. In the follow- 
ing we consider in particular the former case but every- 
thing will hold with appropriate replacements also for a 
measurement of P. 

As shown in figure [S] the set of observables 
2^c"*) Ps"*j PtS can be measured simultaneously by a 
measurement of Stokes component Sy after a 7r/2 rota- 
tion is performed selectively on the sine component of 
the scattered light. The cosine component of the corre- 
sponding photocurrent will give an estimate of and 
the sine component of Multiplying the photocur- 

rent 's sine component by the linear function defining the 
back action mode, equation ljBl|) . will give in addition 
an estimate of Note that the field out of phase is 

conserved in the interaction such that 

„out _ in out _ in ,'-|9^ 

i.e. the results will have shot noise limited variance. It 
is then evident from equation (|5c|l that the respective 
photocurrents together with an a priori knowledge of k 
are sufficient to estimate the mean {X). 

The conditional variances after the indicated measure- 
ments are 

AX2|{,ou.,,„..,^o„.} = (AX-)2-^, (13a) 
Ap2|{,„„._p„..,^o„tj = (AF-)2^^, (13b) 

corresponding to a pure state. Obviously the variance 
in X is squeezed by a factor (1 + k^/2)~^. Note that 
the squeezing achieved in a QND measurement without 
magnetic field but otherwise identical parameters is given 
by {1 + K^)~^ . From this we conclude that the quality of 
the estimate for (X), as measured f.e. by input-output 
coefficients known from the theory of QND measurements 
, can be the same as in the case without Larmor 
precession albeit only for a higher coupling k. 

Equations (|13|l are conveniently derived by 
means of the formalism of correlation matrices 
[3^ . For the operator valued vector R = {X, P, 
Xc,Pc,Xs>PsiXc,i,Pc,i,x^j^,p^^i) equations ©, U^) 
and HB2|) define via _R°"* = S{fi)R"^ a symplectic linear 
transformation S{k). The contributions of ^^'^ 
p™2 to x°^i and as given in (|B2|) are treated as 

noise and do not contribute to the symplectic trans- 
formation S but enter the input-output relation for 
the correlation matrix as an additional noise term as 
follows. The correlation matrix is as usually defined 
by 7i,j = iv{p{RiRj -\- RjRi)}. The initial state is 
then an 10 x 10 identity matrix and the final state is 
^out _ s{k)S{k)'^ + 7noiso where the diagonal matrix 
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B 



EOM x/2 

□ D 



FIG. 5: Scheme for spin measurement: After the scattering a 
7r/2 rotation is performed on the scattered hght modulated at 
the Larmor frequency such as to affect only the sine (cosine) 
component. Standard polarization measurement of Sy and 
appropriate postprocessing allows to read out the mean of 
X (P), leaving the atoms eventually in a spin squeezed state. 



7„„i,^ =diag[0,0,0,0,0,0, 1,0, 1,0](k/2)V15 accounts 
for noise contributions due correlations to second order 
back action modes c.f. equations ljB2|l . In order to 
evaluate the atomic variances after a measurement of 

out out c 
c J -Fs T Ps 



the correlation matrix 7out is split up 



into blocks, 



7out 



A C 
B 



where A is the 2x2 subblock describing atomic variances. 
Novi', the state A' after the measurement can be found 
by evaluating [s^ 



A' = A- lim C- 



1 



T + B 



where T = diag[l/a;, x, x, l/x, n, n, x, l/x] corresponds to 
the measured state. Note that the limit n oo, i.e. 
the projection of the unobserved mode x^2tPc2 onto 
the identity, does not need to be taken explicitly since, 
remarkably, the atomic state after the measurement de- 
couples form this mode. The conditional variances in 
equation H13|l are then just (half the) diagonal entries of 
A'. 



V. CONCLUSIONS 

In conclusion we have presented a simple and realistic 
protocol for telcportation of a coherent state, carried by 
a propagating pulse of light, onto the collective spin of an 
atomic ensemble, a suitable stationary carrier of quantum 
information of continuous variables. The scheme can be 
implemented with state of the art technique and allows 
to surpass any classical strategy for the transmission and 
storage of coherent states under realistic experimental 
conditions. The basic resource in the protocol is a mul- 
timode entangled state as it results form the interaction 
of light with atoms in the presence of an external mag- 
netic field. We showed that Larmor precession enhances 
the creation of entanglement quantitatively and qualita- 
tively. Though the interaction is not of QND type it is 



still possible to perform a state readout on the atomic 
spin as well as to create significant spin squeezing. We 
expect that a proper tailoring of the Larmor rotation 
with time dependent magnetic fields would open up in- 
teresting possibilities to further enhance the creation of 
entanglement and to deliberately shape scattering modes. 

We acknowledge funding from the EU under 
project FP6-511004-COVAQIAL and support from 
Kompetenznetzwerk Quanteninformationsverarbeitung 
der Bayerischen Staatsregierung. 



APPENDIX A: HAMILTONIAN 

In this appendix we present a short derivation of the 
basic Hamiltonian ^ in order to introduce the notation 
used throughout t his pap er. More detailed descriptions 
can be found in [11 El El IHIll Hills IP- The Hamil- 
tonian of the system is given hy H = Hat + Hu + V 

where the atomic part Hat = X^i accounts for the 
external magnetic field along x causing a ground state 
Zeemann splitting of Ml, Hu is the free space Hamilto- 
nian for light and the interaction term V is the level shift 
operator [2^l33| 



V 



(Al) 



which appropriately describes the interaction of off- 
resonant light with atoms. We assume here implicitly 
that the electric field contains only frequency compo- 
nents within a bandwidth h around the carrier frequency 
LOc of the off-resonant coherent probe pulse satisfying 
^ <C Ai?' where the detuning Aj?/ — uj^ — ujf,f'- 
The atomic polarizability density tensor introduced in 
equation IjAlfl is 



a{f) — ^ a*^*'(5(f- 



(A2) 



where r'*^ is the position of atom i. The single atom 

ground state polarizability a consists in general of a 
scalar, vector and tensor part. 



(f (aol + aiFx +a2T 



where d is the relevant reduced dipole matrix element 
of the probed transition, 1 is the 3x3 identity matrix 
and Fx has to be understood to give the vector cross 
product of F with the vector to the right. Each of the 
coefficients aj is a sum of contributions from transitions 
to all excited states manifolds F' . If the detuning is much 
larger than the typical excited states' hyperfine splitting, 
A ^ Whfs, one finds that a2 ^ such that the second 
rank polarizability can be neglected. In this case 



a 



h{A-ir/2) 



(A3) 
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with real dimensionless coefScients aj of order unity and 
r the excited states' decay rate. The non-hermitian part 
of the resulting Hamilton operator describes the effect 
of light absorption and loss of ground state population 
due to depumping in the course of interaction. In the 
following we will focus on the coherent interaction and, 
for the time being, take into account only the hermitian 
component. The effects of light absorption and atomic 
depumping are treated below. 

Coherent interaction Since scattering of light occurs 
predominantly in the forward direction 221 it is legiti- 
mate to adopt a one dimensional model such that the 
(negative frequency component of the) electric field prop- 
agating along z is given by 

E^-\z,t) = E'^-\z)ey+£^-\z,t)e^ 
E^-^z) = p{u;c) [duja\uj)e-''''' 

Jb 

£^-\z,t) = pK)y^4^7v7re-'('=^^-"^*) 

where p{uj) = ^ hwd ^-i^tQAc and A denotes the pulse's 
cross sectional area, the overall number of photons 
in the pulse and T its duration. We restrict the field 
in X polarization to the classical probe pulse, since only 
the coupling of atoms to the y polarization is enhanced 
by the coherent probe. Furthermore we implicitly assume 
for the classical pulse a slowly varying envelope such that 
it arrives at 2; = at i = and is then constant for 
a time T . Combining this expression for the field with 
expressions (|A2|) and ljA3|l for the atomic polarizability in 
equation IjAip yields 

VAttJT Jb J ^ 

where we defined a dimensionless coupling constant 
K = \J NphJuJcai(P / hegcAA and an atomic spin density 

jz{z) — J2i Fz^^S{z — z^')). In this expression we skipped 
terms proportional to oq which will give rise only to a 
global phase shift and included for convenience the square 
root factor with J = NatF where Nat is the number 
of atoms. Note that the coupling can be expressed as 
K = \/NphJaiaT /2A/S. with a the scattering cross sec- 
tion on resonance. 

We now define field quadratures for spatially localized 
modes 0,0 as 

xiz) = (dojiaiuj^e-'^^"-^^'- +h.c\ , (A4a) 

\/47r Jb ^ ' 

p{z) = --^ Jduj (a{Lj)e-'^''^-''^'' - h.c)j {A4h) 

with commutation relations [x{z),p{z')] — icS{z — z') 
where the delta function has to be understood to have a 
width on the order of c/b. Since we assumed that ft b, 
the time it takes for such a fraction of the pulse to cross 
the ensemble is much smaller than the Larmor period 



During the interaction with one of these spatially 
localized modes the atomic state does not change appre- 
ciable and we can simplify the interaction operator to 

V = hK{JT)'^/^Jzp{Q) where = Y^i^z^ and we as- 
sumed that the ensemble is located at z = and changed 
to a frame rotating at the carrier frequency ojc- 

A last approximation concerns the description of the 
atomic spin state. Initially the sample is prepared in a 
coherent spin state with maximal polarization along x, 
i.e. in the eigenstate of Jx with maximal eigenvalue J. 
We can thus make use of the Holstein-Primakoff approx- 
imation [2^ which allows to describe the spin state as 
a Gaussian state of a single harmonic oscillator. The 
first step is to express collective step up/down operators 
(along x), J± = Jy ± iJz, in terms of bosonic creation 
and annihilation operators, [b, 6^] = 1, as 

J+ = V2J^1 - 6t6/2J6, J_ = V2Jb^^t - 6t6/2J. 

It is easily checked that these operators satisfy the 
correct commutation relations [J_|_,J_] = 2Jx if 
one identifies = J — Wb. The fully polarized 
initial state thus corresponds to the ground state 
of the harmonic oscillator. Note that this map- 
ping is exact. Under the condition that {b%) <C J 
one can approximate J+ ~ \/2Jb, J-_ ~ •\/2J6^ and 
therefore Jz — —i\/ J /2{b — b^). Introducing atomic 
quadratures X = {b + b^)/V^ and P = -i{b - h^)/^/2 fi- 
nally yields the desired expression for the interaction 

V = hnT-^/^PpiO). 

In terms of atomic quadratures the free Hamilto- 
nian for atoms is H = h£l/2{X'^ + P^). In the 
frame rotating at the carrier frequency the action 
of Hii on the light quadratures x{z),p{z) is simply 
i/h[Hii,x{z)\ — —cdzx{z) and likewise forp(z). 

Noise effects The antihermitian part of the level shift 
operator IjAljl describes depumping of ground state pop- 
ulation and photon absorption. The effect of the lat- 
ter process can - as far as it concerns the performance 
of the teleportation protocol - be treated on equal foot- 
ing with mode mismatch and finite detector efficiency. 
This is done in section IIII Bl Loss of ground state pop- 
ulation on the other hand will eventually cause degrad- 
ing of atomic polarization due to spontaneous emission 
events. For a single atom the dominating term de- 
scribing this process stems from the scalar part of the 
polarizability and is given by Vioss = ihrj^/AT where 
'q = NphaoWcTd? /2hI^^eoAc = NphaoaT'^ lAA/:^'^. It is 
possible to have 77 <C 1 and at the same time a large 
coupling K ~ 1. For a thermal cloud of atoms an addi- 
tional source of decoherence are light assisted collisions 
which in fact dominate the decay process. Assuming a 
transverse relaxation at an overall rate (3/T with (3^-2 
the exponential decay during the interaction can to a 
good approximation be treated linearly which leads to 
equations H1U|I . The Langevin noise operators fx.p can 
in prin ciple be derived by a microscopic model as is done 
in 1/1 for dephasing due spontaneous emission. 
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APPENDIX B: BACK ACTION AND INPUT 
MODES 

Back action We evaluate here the input /output rela- 
tions H5cl5d|) for the cosine and sine components of the 
in-phase field. For the former we take equation (j3d|) at 
^ — CT, t = T, multiply by •y/2/T cos(r2T) and integrate 
over T from to T. Using equations H3al I3b|l and the 
approximate orthogonality of cos{^It) and sin(riT) one 
finds 



V2 



,2 , dr / dT'[cos(17r)2sin(f7T')p(cr',0) - 
T ' Jq Jo 

— cos(17t) sin(riT) cos{V,t)p{ct' , 0)]. 

After interchanging the order of integration, 
Jp^dr /Jdr' — > /q^cIt' /^dr one can perform the in- 
tegration over r. Neglecting all terms of order n^^ or 
less where uq = flT 3> 1 one finds 



4^' 
V2 



r3/2 



■ sm{flT)p{cT, 0). 



The last term represents back action of light onto itself. 
It can be expressed as a sum of two terms, one propor- 
tional to p™ and another one proportional to 



Ps,l 



3/2 



T 



dr ( - - r ) sin(17r)p(cr, 0). (Bl) 



It is easily verified that the back action mode defined 
by this equation and the corresponding expression for 
2. back jg canonical [a;s"i,Ps"i] — «[1 — ©(rig ^)] — « and in- 
dependent from all the other modes introduced so far, 
f.e. [a;™,Ps"i] = ©(tiq^) ~ 0. The variance is thus 
(Ap™-^)^ — 1/2. Repeating the calculation for a;°"' with 
appropriate replacements and a definition of p J,"]^ analo- 
gous to equation (|Bip finally yields equations H5cl5d|l . 

In a similar way input-output relations for the back 
action modes itself are derived. In particular for the in 
phase components one finds 



1 / K 



Pc 



15 V2J (^2^) 



PTa- (B2b) 



15 V2 



In both equations the third terms on the right hand side 
describe contributions of second order back action modes 
defined by 



Ps,2 



-1^ 



1/2 i-T 



drl Tt 

6 



sin(r2r)p(cr, 0) 



and similarly for x™2 and the cosine component. These 
modes are again canonical and independent. As a side- 
mark we note that, formally, it is possible to define scat- 
tering modes of arbitrary order whose mode functions are 



given in general by products of Legendre polynomials and 
cos{Q,t) [sin(rit)] resulting in a hierarchy of input-output 
relations similar to l|B2(l . 

Input state The input field, propagating along the 
positive y direction and polarized along z (see figure |21l, 
is described by operators \b{ijj),h^ [u')] — 5[ijJ~uj') in fre- 
quency space and [y{y), q_(y')] — ic6{y — y') in real space. 
{y is the quadrature operator for the field in-phase and y 
on the other hand is the position along the y-direction.) 
In analogy to equation (|A4I) these bases are connected 
via 



1 



47r J6 
i 



duj (b{uj)e-'^''''-'''>y + h.c^ 
' d.(6He-(^^-)^-/..c, 



viv) = 
q{y) = 

As shown in section UlI Al we can teleport the mode 



y 



9c 



Qsin) 



and the claim is that this corresponds to a pulse at the 
upper sideband ujc + ^- Using the definitions of cosine 
and sine modes (0J as well as y(cr, 0) = y(— cr, 0) and 
the same for q{y) we have 

y = —j= I dr [sin(r2r)?;(cT, o) -I- cos(r2T)5(cr, 0)] 
V T Jo 

= f dr fduj\b(uj)e'^'^^+^^-'^'>^ - h.c. 

Jo Jb ^ 

q = ^ / dr [ cos(f2T)j/(cT, o) — sin(f7r)g(cT, 0)] 

6(tj)e*('^'^+^^-")^ + h.c. 



-1 



^ f 
dr duj 

Jb 



To explicitly see that this corresponds to a pulse centered 
at the upper sideband it is convenient to change to a 
more precise model by replacing the I/VT factor, which 
is just the pulse's slowly varying amplitude function in 
a simple square well approximation, by a function A{t) 

of dimension s"^/^ normalized such that jQ'dr|A(r)p — 

1. Its Fourier transform A{lu) — --^ jjdryl(r) exp(iu;r) 

is centered at zero and has a width l/T = Alu ^ 51 
in accord with our condition 1 <JC flT. Replacing now 
1/Vt by A{t) (of course inside the integral over r) in 
the expressions for y and q yields 



6(w)e'('^'^+"-")^ - h.c. 



An 
i 



dryl(r) dcj 



^/2 

-1 



duj 



Lu)b{uj) — h.c. 



q = —= duj A{ujc + — uj)b{uj) + h.c. 
v2 Jb 



This is evidently a mode whose spectral mode function is 
the same as the classical pulse but is centered at ujc + ^. 
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APPENDIX C: FEEDBACK 

The feedback in continuous variable quantum telepor- 
tation is sometimes described by equations equivalent 
to ||SJ) but with a classical random variable describing 
the measurement outcome in place of the operators cor- 
responding to the chosen displacement, which - though 
giving the right result - is mathematically questionable. 
We point out that relations © hold stricto sensu as op- 
erator identities. This is true for mixed and even for 
non-gaussian states, as we will show below. 

Consider a bipartite system of A^ + l modes and denote 
the first mode as system A and the remaining N modes 
as system B. Let the state of the compound system be 
given by pab- Our aim here is to describe protocols 
which consist of the following steps: 

Measurement On system B a set of commuting ob- 
servables {fi, . . . ,r]y} is measured where each of the op- 
erators fi is either Xi or pi, one of the quadratures of 
mode i in B. Let the corresponding measurement out- 
comes ri be arranged in a vector = (ri, . . . , tat). 
With the eigenvalue equation fi|ri)s = ri\ri)B, where 
\ri) B is the generalized eigenstate of fi, the normalized 
state of system A conditioned on the measurement out- 
comes is 

Pa\^) = Bin, ■ ■ ■,rN\pAB\ri, . . .,rN)B/p{R)- 

p{R) is the probability to get the measurement outcomes 
R and is normalized as Jd^rp{R) = 1. 

Feedback Depending on the measurement outcomes 
system A is then displaced in xa and pA by an amount 
Rg^ and Rg^ respectively where gx{p) are any real N 
dimensional (row) vectors determining the strength with 
which each outcome is fed back into system A. In tele- 
portation protocols these coefficients are usually referred 
to as gains. The state of system A after the feedback 
operation is then 

pf{R)^D\p^]^\R)DA. 

Da = DAiRg^, Rg^) = eyip{iRg'^pA) e^p{-iRg'^XA) is 
the unitary displacement operator effecting the de- 
sired transformations DaxaDI^ — xa + Rg^ and 

DAPADl=pA + Rg^. 

Ensemble average On average over all measurement 
outcomes, weighted with their respective probabilities, 
the state of system A is 

PA= jd''rp{R)pf{R). 



Combining this with the expressions for and pj^ 
above one can express 



PA = yd rD\B{ri,-..,rN\pAB\ri,.-.,rN)BDA 
= Jd^rsiri, . . ■ ,rN\Dl^gPABDAB\ri, ■ ■ ■ ,rN)B 

= tTBlD^gPABDAB} (CI) 

where the trace in the last line is now taken with respect 
to both system A and B. In going from the first line to 
the second line we made use of the identity 

\ri, ■ ■ . ,rM)BDA = DabVi, ■ ■ . ,rAr)s (C2) 

with the unitary operator Dab defined as Dab — 

DAB{Rg^,Rg^) = exp{iRg^pA)exp{~iRg'^pA) where 

R — (fi, . . . ,fjv) is now the vector of operators ff and 
Dab acts on both systems A and B. Note that identity 
(|C2|) is valid only for commuting observables ff . The 
resulting equation ljCl|l is the key point in this consider- 
ation. 

Observables in the ensemble average Consider finally 
the mean of, for example, xa after the measurement and 
feedback procedure, i.e. with respect to the ensemble 
averaged state pA- It is given by 

{xa) = trAixApA) = ^T:AB{DABXAD\gPAB} 

= trAB{{xA + Rg^)pAB}- 

From this identity and the corresponding expression for 
(pa), which both are true for all initial states pab, we 
can deduce the operator identities 

xi" = XA + Rgl , p^A =VA + Rgp 

where x^,p^ describe the final state of system A in 
the Heisenberg picture and means have to be taken with 
respect to the unchanged initial state of both systems 
A and B. If pab is a pure Gaussian state the last two 
equations fully determine the final state pA- This was 
used in equation (jS)). Note that these considerations are 
easily extended to situations in which system A consists 
of more than one mode. 
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